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ABSTRACT 



The solar atmosphere is structured and inhomogeneous both horizontally and 
vertically. The omnipresence of coronal magnetic loops implies gradients of the 
equilibrium plasma quantities like the density, magnetic field and temperature. 
These gradients are responsible for the excitation of drift waves that grow both 
within the two-component fluid description (in the presence of collisions and 
without it) and within the two-component kinetic descriptions (due to purely 
kinetic effects). In the present work the effects of the density gradient in the 
direction perpendicular to the magnetic field vector are investigated within the 
kinetic theory, in both electrostatic and electromagnetic regimes. The electro- 
magnetic regime implies the coupling of the gradient-driven drift wave with the 
Alfven wave. The growth rates for the two cases are calculated and compared. 
It is found that, in general, the electrostatic regime is characterized by stronger 
growth rates, as compared with the electromagnetic perturbations. Also dis- 
cussed is the stochastic heating associated with the drift wave. The released 
amount of energy density due to this heating should be more dependent on the 
magnitude of the background magnetic field than on the coupling of the drift and 
Alfven waves. The stochastic heating is expected to be much higher in regions 
with a stronger magnetic field. On the whole, the energy release rate caused 
by the stochastic heating can be several orders of magnitude above the value 
presently accepted as necessary for a sustainable coronal heating. The verti- 
cal stratification and the very long wavelengths along the magnetic loops imply 
that a drift-Alfven wave, propagating as a twisted structure along the loop, in 
fact occupies regions with different plasma-/? and, therefore, may have differ- 
ent (electromagnetic-electrostatic) properties, resulting in different heating rates 
within just one or two wavelengths. 
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Introduction 



Observations and theoretical studies in the past 70 years have dramatically increased 
our knowledge and understanding of the physical processes in the solar atmosphere. 
However, the basic starting puzzle of the problem of coronal heating still remains elusive in 
spite of the obvious progress made in the domain. In fact, new data collected in the course 
of decades have additionally increased the complexity of the problem as more and more 
fine details related to the heating have emerged. These include the preferential heating of 
plasma particles in the dir ection perpendicular to the magn etic field vector, resulting in 



a temperature anisotropy (JLi et al 



1998 : 



Cuseri et al. 



19991 ). and the preferential heating 



of heavier particles. As a matter of fact, heavier ions appear to be hot ter than lighter 



ions, the latter on the other hand appear to be hotter than electrons (jCranmer et al. 



20081). Moreover, extre mely strong electric fields (above 100 kV/m) have been detected 



( IZhang fc Smartt 



19861 ). Those electric fields accelerate particles and, in general, the 



distribution functions of the plasma species in the outer solar atmosphere ca n considerably 



be different from a Maxwellian distribution ( iVasyliunas 



1968 



Cranmer 



19981). 



In our recent papers (iVranjes fc Poedts 



2009a 



d|) a novel approach and a new 



paradigm for the coronal heating has been put forward. The model is based on the drift 
wave theory, a well known subject in the general pl asma t h eory, in laboratory plasma 



KellevI (119891). Yet, the drift wave 



physics, and even in terrestrial ionospheric research 
theory is completely overlooked in the context of solar plasmas. It implies the abundance 
of free energy for the instability of the drift wave already in the corona. That energy 
is stored in the gradients of the density, temperature and magnetic field. It makes the 
waves growing and it results in heating due to the polarization drift effects. The heating 
is stochastic by nature; for short enough wave-lengths in the direction perpendicular to 
the magnetic field vector, gyrating plasma particles feel a space-time variation of the 
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wave-electric field, and their motion becomes stochastic and equivalent to the increase of the 
temperature. The nature of the heating is such that it acts essentially in the perpendicular 
direction, and, in addition to this, more massive partic les are in fact mo r e effect ively 



Vranies fc PoedtsI (J2009a|jb 



Idjdl) 



heated by that mechanism. The analysis performed in 

was focused on the electrostatic domain of the drift wave instability. This implies a small 

plasma-/? = 2fiQnQK,T/BQ, e.g. of the order or below me/mi. Howev er, even in that domain 



the plasma may support electromagnetic (EM) perturbations too (JKrall 



19681), although 



as a rule those will not be well coupled to the electrostatic ones. For a plasma-/3 value 
above rrig/wij, the coupling will effectively take place and, as a first manifestation of this, 
only the bending of the magnetic field vector may be taken int o account. Such a couple d 



Vranies fc PoedtsI fl2006h 



drift-Alfven mode has in fact been studied in our earlier work 
by using two-component fluid theory, with a complete and self-consistent contribution of 
hot ions effects, appropriate for the hot solar corona. Within such a two-fluid theory, the 
drift-Alfven mode is destabilized in the presence of collisions and, for a large enough parallel 
wave- number, the mode has all the features of a growin g Alfven wave. This is because o ' 



an ex change of identities of the drift and Alfven modes (jWeiland 



2000 



Vranjes &: Poedts 



20061 ) occurring in a certain parameter domain. 



On the other hand, the drift-Alfven wave instability within the collision-less kinetic 
theory has a completely different nature, and this will be the subject of the present work. In 
particular, we shall investigate the difference in the growth-rates of the electrostatic (drift) 
and the electromagnetic (drift-Alfven) modes. Such a difference is expected because of the 
following two opposite effects: i) a part of the energy that drives the instability is spent 
on the bending of the magnetic field vector and this should in principle reduce the growth 
rate, yet in the same time, ii) this bending should partly reduce the electron mobility in the 



parallel di rection. The effect s of su ch a reduction should be similar to electron collisions 

2~ 



studied by 



Vranjes &: PoedtsI ( 120061 ) and, as a result, the growth-rate may become increased. 
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Hence, the total outcome of the EM effects will then depend on the mutual ratio of these 
two opposite effects. A local analysis will be used. This is very appropriate for a geometry 
without magnetic shear, and for the case in which the perpendicular component of the 



wave- 



f jKrall 



ength is much shorter than the characteristic lengths of the equilibrium gradients 



19681) 



2. Basic equations 

In the case of a plasma beta in the range rrie/mi < /3 ^ 1, it is appropriate to take 
into account only the bending of the magnetic field . The perturbed density for the species 



j = e,i is, in that case, described by (IWeiland 

qjUjo 



2000): 



Ujl 



^mih) 



kTj 



X 



W 



1 + {u - ui)(f)iy^ — - 

U2 — rriQj 



U2 — mVtj 



I i^Z I '^Tj 



- 1 



+(^i-0i)(l-^) Y,K^{hj)W 



UJ2 — rriQj 
oj / '' — ' ~ "' \ \kAvTj 

Here, cui = uj^j — kyg/Vtj, uj2 = uj + kyg/Vtj, uo^j = v^jky, where w^j = — e^ x Wpjo/igjUjoBo) 

is the diamagnetic velocity. The equilibrium magnetic field and density gradient are 

Bq = BqCz and Vn^o = —exdrijo/dx, and we use a local approximation and Fourier 

analysis with small perturbations of the form ~ f{x)exp{—iujt + ikyy + ik^z), where 

f{x) is the x-dependent amplitude and \d/dx\ -^ \ky\. In the terms ui,2 we have 

the gravity effects that are included through the Maxwellian distribution function 

with a rrijgx term in the exponent. Below, it will be kept for ions only. The other 

notation is as follows bj = kypj, Am{bj) = Im{bj)exp{—bj), pj = Vtj/^j, v'^j = nTj/nij, 

flj = qjBQ/rrij, Im is the modified Bessel function of the first kind and order m, 

^(x) = (27r)-'^/^ J_ °°?7exp(— ?7^/2)c??7/(r7 — %). The terms and i}) describe the potential of 
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the electric field (IWeiland 



20001), E = —W±(j) ~ Czd^jj/dz, and for that reason an additional 
equation for the parallel current is needed in order to have a closed set 

'u2 — mftj 



Jjzl 



q\o_ 



{U - UJi)<Pi^Am{hj)W I 



\ I "^2 I '^Tj 



+(V^i - 0i) (l - ^) X 



X 



^K^{hj){uj2~rnnj)W 



UJ2 — mQj 

I f^Z I "^Tj 



(2) 



For electrons it is good enough to use a negligible mass limit so that Ao(&e) — 1? while 
the deviation from unity of the corresponding term for ions is a finite Larmor radius 



ow fij, one keeps only t 



le term m 



1992 



Weiland 



2000 : 



Bellan 



20061). 



effect. In the limit of frequencies 10721 much be 

in the summation for both electrons and ions ( iStix 

For similar reasons, in the case |x| < 1, we shall use the approximate expression (for 

electrons) W{x) — 'i(7r/2)-'^/^x6xp(— x^/2) + 1 — x^ + X'^/4 . . ., and for ions |x| > 1, 







w{x) ^ lijilM^ 



becomes ( IWeiland 



yexp 



20001 ) 



-X^/2) — 1/x^ — 3/x^ + . . .. For these two species Eq. ([T]) then 



no 



nTf. 



(Pi + ^(Pi\-) TTI 6xp 



UO 



2klvl, 



+ (^i-0i)(l-^)[l + ^( 



'7r\i/2 



X- 



UJ 



\i^z\VTe 



exp 



uo 



2klu 



2 

z ^Te 



(3) 



no 



e 



01 -01 ■ — Ao(Oi 



7r\i/2 UJ2 



\kAVrj 



002 
exp 



1 _|_ ^Z^Tt _|_ ^^Z^Ti 



ujI 



^2 



Un 



2klvl, 



+ (V^i -0i)x 



X (1 - ^) A„(. 



7r\i/2 UJ2 



exp 



wl 



ZrC^Vrj^. 



, ,2 






(4) 



3. Electrostatic drift wave instability 

In the electrostatic limit, one may set 0i = ipi in the above expressions. Using the 
quasi-neutrality condition rii = rif. and Eqs. fl3]H|) . one directly obtains the dispersion 
equation in the form i?eA(ci;, k) +iI'mA{uj, k) = 0. The frequency is assumed to be complex, 
in the form u = Ur + i'y. Setting ReA{ujr, k) = 0, one then obtains the following equation 
for the real part of the frequency: 

-Ao{bi)klclu}rU}l^ + Ao{bi)uiklclul^ - 3Ao{bi)kjvliclur 



+ 3Ao(6i)A;XiC>i = 0. 



(5) 



In the limit of a negligible ion response along the magnetic field vector, \ur/kz\ ^ c^ with 



n Tp/mj, and for sma l 



used in 



Vranjes fc PoedtsI (120093 



gravit y effects, this gives the electrostatic drift wave frequency 



mm 



UJr 



uJ*iAo{bi 



W*eAo(&i 



Here, w* 



1 - Ao(6.) + r,/T, 1 + [1 - AomrjT, 

-u^eTi/Te. Observe that after setting Ao(6j) — 1 — &«, the term in the 



(6) 



denominator becomes 1 + fc,jp:, where 



comparison and as a guideline (jBellan 



2006 



c.s/O,;. Using the two-fiu id description for 



Vranjes fc Poedtsi 120061 ). it can be shown 



that this same expression (describing the finite ion inertia) directly follows from the ion 
polarization drift, the latter playing an esse ntial role in the process of stochastic heating 



( JMcChesney et al. 



1987 



Sanders et al. 



19981 ) that will be discussed below. 



On the other hand, for a purely parallel propagation from Eq. ([5]) one obtains the 
ion-acoustic (lA) mode in plasmas with hot ions 



^r - klcluj'^r - ^ktvliCl 



0. 



(7) 
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The solutions cof. = {k'ici/2)[l + (1 + 12Ti/TeY/^] for Ti = T^ yield the frequency 

Ur — ±1.5 fczCs- Hence, the ion temperature plays no important role in the real part of the 

I A wave frequency, while the opposite is certainly true with its imaginary part. 



The growth rate 7 is obtained approximately from 7 
This yields 



-ImA{uJr,k)/[dReA/dLo] 



7 = - 



u. 



2r 



77 \ 1/2 

2/ a;*eAo(&j 



UJr - ^* 



+Ao(Oi) — -^7-77 — exp 



/!r2 2 

1 + ^(3 

W2r 






exp 



^Ir 



w: 



2klvl, 



2k^Vj,^ 






Au' 



'2r 



The frequency on the right-hand side in Eq. (|8]) is to be obtained from Eq. (|5]). Note that 
in the absence of gravity uj2r -^ ^r-, and in addition, for a negligible ion response along t he 



magnetic field vector, Eq. ^ yields the growth rates from IVranjes fc PoedtsI ( l2009al JblJdJdl). 
For the assumed geometry u^i is negative and the ion contribution will always tend to 
reduce the growth rate, while electrons will make the mode growing. Clearly this purely 
kinetic instability can only take place provided that the frequency Ur is below uj-^e- 



4. Electromagnetic perturbations 



To proceed with the electromagnetic perturbations, a procedure similar to the 
derivation of Eqs. (E]) and (jl]) yields the electron parallel current 



]ezl 



T, 



e^riQ 



kzKl e 

7r\i/2 OJ 



{U -U^e) (1 + ^Te) ^"1, 



(9) 



. Z / I fvj; I fye 

The corresponding expressions for the ions are: 



exp 



u 



2klvl, 



Jizl 






[u - uji)Ao{bi] 



1 + — U-i ^01 

00 / 00 . 
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X 



1 H — I - iTi 



I [2 
UJ2 



UJ2 



(10) 



T, 



7r\l/2 UJ2 

2/ |A;r|fT 



exp 



w; 



2klvl. 



w« 



M 
f^r 



The first necessary equation is obtained as above, by using the quasi-neutrahty condition 
ni\ = riei and Eqs. ([3]and|l]). The second equation follows from the Ampere law that, with 
the help of Eqs. (ED and (flO]), yields 

Sl + fl 



Ipl = 01 



S2 + /2 



;ii) 



si = kyp^k.c^, S2 = si + 



UJ[UJ*e 



c^)(l + ^T, 



z 



UJ2 



/i = —ko{hi)ujg{uj - ui) (a - iTi) , /s = /i — 



Ua 



a 



ul 



,2„,2 



'^z'^Ti / 1 _|_ '^z'^Ti 



ul 



Here, c^ = Bo/{fionQmi). The terms /i,2 originate from the ion parallel current and in many 
situations can be neglected. 



The potential (jTTj) is used to eliminate ipi in the quasi-neutrality condition and in the 
end one obtains the following dispersion equation: 



'T,; 



UI2 



Ao{h){l + a-iTi 



+p ^^-^^ - Ij —^Uh) (.T„ - a 



:i2i 



Here, 



T. 



TT 



1/2 



W — W* 



exp 



u 



2kH 



2 

2 '-'Te 



The parameter p is set here by hand, and for convenience only; taking p = is equivalent 
to the electrostatic limit discussed in the previous section, while p = 1 implies the 
electromagnetic perturbations that are of interest here. 
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The dispersion equation flT2|) describes the coupled Alfven and drift waves as well as 
the ion parallel (acoustic) response, together with the gravity and finite Larmor radius 
effects for ions. Note that, for the parameters used further in the text, the gravity drift 
frequency kyg/Qi is usually negligible. 

The wave spectra are obtained following the same procedure as before yielding from 
the real part of Eq. flT2|) : 



= 1 + 



T, 



l-^^^^-^^Ao(&.)(l + «. 



UJ2r 



. -Sir + fir _ ^ 



^. ^1 1- — -Ao(&.K7^ 1-- 

,S2r -r J2r / [ ^r -'i \ ^'^ 

Here, the index r denotes the real part of the corresponding expressions. 



w* 



T. 



Ui 



(13) 



It is interesting to compare these derivations with the results from the two-component 
fiuid theory. Omitting the ion parallel response and gravity, Eq. (IT3l) yields 

(14) 



(w - w,e) (w^ - ^*i^ - klcl) = klclk^pliu - u^i). 



This equation is 



e xactly the same as the corresponding two-fiuid equation from 



Weiland 



Vranies &: Poedtd (120061 ) . and it is also obtained from the kinetic derivation in 
( I2OOOI ). We stress that such a perfect agreement between the two (fiuid and kinetic) 
descriptions is only possible if the two-fiuid derivations self-consistently include the 



gyro-v i scosity stress tensor contributio ns. Details on these issues can be found in 



Weiland 



fcoooh 



Vranies fc PoedtsI ( 120061 



2009el ). From Eq. (TT4l) it is seen that the coupling between 
the drift mode u = u-^e and the Alfven mode is due to the right-hand side, which here 
appears due to the finite- ion-mass effect k'lpl. In the fluid description, the latter term 
originates from the ion polarization drift Vip = — (c}/9t)[V^0i/(f2ji?o)]. 
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5. Growth rates 

In both instabihties discussed in the previous section, the perpendicular ion motion 
is essentially the same: the typical dominant velocity is due to the E x i?-drift and this 
holds as long as Xy ^ pi. The latter condition may be relaxed, the relative contribution 
of the polarization drift in that case increases and the foreseen stochastic heating for the 
given electrostatic-drift and drift- Alfven in stabilities will hav e a sim ilar nature. Note 



that experimental verification performed in 



McChesney et al. 



(119911 ) in fact involved the 



excitation of the drift-Alfven waves. It may be of importance to check the growth rates of 
these instabilities, calculated for the same or a similar set of physical parameters. This may 
give the answer about their relative importance in the solar corona. 

To check our model and the differences between the electrost atic (ES) and 



electrom agnetic (EM) cases, we take a set of parameters similar to 



Vranjes fc Poedts 



(12009 tf ): no = 10^^ m^^, Bq = 10"^ T, and T^ = Ti = 10^ K. We further set 
Ln = [{dnQ/dx)/nQ]~^ = s ■ 10^ m, and take the parallel wave-length A^ = s ■ 10^ m. 
The plasma (3 for the present case is 0.64me/mi; as shown below this can be taken as an 
appropriate electrostatic domain. The parameter s can in principle have any value (e.g. in 
the interval 1 — 10^), only bearing in mind the necessity of staying reasonably well within 
the previously imposed conditions used in the derivations. The simultaneous variation 
of Ln and A^ by changing the factor s is introduced for convenience only: as shown in 



Vranjes fc PoedtsI (l2009bl ) for the electrostatic limit, by doing this it turns out that the 
ratio 'j/ur remains exactly the same regardless of the value of s, while the actual values of 
both quantities 7 and Ur can, in fact, drastically change. 

In Fig. 1 we plot the frequency Ur for the case s = 1, and in terms of the perpendicular 
wavelength Xy. The four lines in the figure represent the real part of the frequency calculated 
from Eq. (iT2l) in the following manner. The full line is obtained after setting p = 0, that is 
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equivalent to simply neglectin g the EM effects. I n that sense it corresponds to the purely 



electrostatic analysis given in 



Vranjes fc Poedtj ( l2009bl ). The dotted line is obtained after 



setting p = 1, yielding consequently the drift wave frequency when the EM effects are taken 
into account. As expected, in view of the given small plasma-/3, the frequency remains 
almost unchanged. The frequency is passing through a maximum, and this follows from the 
fact that Ur ~ ky/{l + kypl) (note that kyP^ ~ 9 at A^, = 0.2 and kyP^ ~ 0.36 at A^, = 1). 
The plasma-/3 can be varied by changing several parameters. In the present case this is 
done by the variation of the number density. Hence, we take it to be wq = 2 ■ 10^^ m~^ and 
uq = 5 ■ 10^^ m~^, and calculate the frequency from Eq. flT^ . This is represented by the 
dashed and dash-dotted lines, respectively. It is seen that the drift-wave frequency becomes 
reduced and this can be attributed to its coupling with the Alfven wave. Note that setting 
s = 10^ (thus simultaneously changing L„ to 10^ km, and A^ to 10^ km, implying larger 
coronal loops) reduces the frequency approximately to Ur/s, while in the same time the 
dispersion lines keep exac tly the same shape. Hence, similar to the electrostatic analysis in 



Vranjes fc Poedtd (J2009bl ). the reduction of frequency due to variation of s remains more or 



less the same even in the presence of EM perturbations discussed here. 

The plot of the growth rate corresponding to the frequencies from Fig. 1, is given in 
Fig. 2. It shows that the EM effects can make an important modification of the drift wave 
for larger values of the plasma-/3. The reason for the reduced growth rate for short Xy 
can be understood from Eq. (IT^ : for larger ky the coupling term on the right-hand side 
becomes more important, a greater amount of energy is spent on the Alfven mode and, 
because of the fixed amount of the free energy stored in the background density gradient, 
th e growth rate is th e refore reduced. Hence, the physics of the stochastic heating predicted 



m 



Vranjes fc Poedtd (l2009al JblJdJcll) will remain nearly the same even in the case of the 
coupling with the Alfven mode, provided a low plasma-/?. On the other hand, for a larger 
plasma-/3 the increasing electromagnetic effects will impose longer growth times. 
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Next, we check the drift mode behavior in terms of the parallel wavelength \z and this 
for several different values of the plasma-/?. For that purpose Eq. flT^ is solved numerically 
and the results are presented in Fig. 3 for the plasma number densities uq = 10^^ m~^, 
no = 6 ■ 10^^ m~^, and no = 10^^ m~^ that yield /3 = O.GAme/mi, 3.8me/mi, and GAme/rrii, 
respectively. Other parameter values are A„ = 0.5 m, L„ = 1 km, Tp = Tj = 10^ K. Th e 



given s hape of the 'j/cOr lines are equivalent to those from Fig. 3 in IVranjes fc Poedts 



(l2009bl ). Here too, the increased EM effects reduce the growth rate. On the other hand. 



similar to IVranjes &: PoedtsI (l2009bl ). for relatively short parallel wavelength components 
the instability vanishes. This is due to mobile electrons which now have to move within 
shorter distances in the parallel direction in order to short-circuit the potential buildup 
caused by the wave. 

The plasma-/? can change also by varying the temperature and/or magnetic field. 
However, the drift wave frequency is proportional to the temperature and also strongly 
depends on the magnetic field [see Eq. ([6])], so that the effect of the perturbed magnetic 
field alone on the drift-wave in that case is not so transparent. In Figs. 4 and 5 we give 
the drift wave frequency and growth rate in terms of the plasma temperature for several 
values of the plasma density. Other parameter values are \y = 0.5 m, A^ = 200 km, and 
L„ = 1 km. The electromagnetic effects are again most effectively seen by taking several 
values of the plasma density and varying the temperature. At T = 1.4 ■ 10^ K the frequency 
is reduced by factor 1.4 for the density increased from uq = 10^^ m~'^ to no = 10^® m~^. At 
the same time the corresponding growth rate from Fig. 5 is reduced by a factor 3.2. 

The graphs presented in this section are solely for the drift wave part of the spectrum 
from Eq. (TT2|) . The Alfven mode, that is also described by Eq. (fT2l) . plays no important 
role in the present study dealing with the stochasti c heat ing. The coupling of the two 



modes is in fact given in detail in 



Vranjes &: PoedtsI (120061 ) using the two-component fluid 
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descriptions. 

In view of the parameters used in this section, it is seen that the quasi-neutrahty 
condition used in the derivations is well satisfied. This because the Debye length A^ is 
typically around 1 mm so that k'^X'^ — 0.0004 ^ 1, and in the same time the equivalent 
condition Q'^/u^^ ~ 0.0005 ^ 1 is also satisfied. For these parameters, the coronal plasma 
from our examples is rather similar to the tokamak plasma. 

6. Application to heating 

Considering a specific single coronal loop and in view of the given geometry that 
implies a very elongated wave front in the axial direction, A^ ^ A_l, the results presented 
above may imply the following. The waves propagate both axially and poloidally, with 
drastically different wavelength components in the two directions. In the cylindric geometry 
of a magnetic loop, the wave front is thus twisted along the loop. An extremely large axial 
component of the wavelength implies a wave that simultaneously takes place in areas with 
gradually varying (with altitude) plasma parameters and consequently different plasma-/?. 
At higher altitudes with a lower density, the perturbations may be electrostatic and develop 
on a shorter time scale. The associated heating may rapidly develop at such places and 
it can then spread along the common wave-front towards lower regions where, due to the 
increased plasma-/?, it is additionally accompanied by the EM effects that develop on longer 
characteristic times. 

However, this scenario with an increased plasma-/? due to higher plasma density can be 
partly counteracted by the lower temperature at lower altitudes, and as a result the energy 
release rate and the heating, together with the variation of magnetic topology, may not be 
so drastically different along the given magnetic loop. One possible example of this can be 
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seen from Fig. 5: for the given normalized temperature T = 1 from the hne 1 (point A in 
Fig. 5), the plasma- /3 is O.GAmg/mi, while for example at T = 0.4 from the line 3 (point B), 
the plasma-/3 is 2.5'm(./mi. This is of course just a rough estimate because the points A 
and B belong to two separate dispersion lines. One particular plasma mode is determined 
by one particular dispersion line, yet in view of the parameters changing with the altitude 
such a transition may be expected. Hence, the drift wave spreading along the loop will have 
an electrostatic nature at the first point, and it will be EM at the second one, implying a 
difference in the heating rate and the magnetic variation. 

The magnetic field intensity too varies with the altitude, and this may additionally 
change both the plasma-/? and also the wave properties. Hence, assuming the magnetic 
field is stronger by a factor 3 (i.e. taking Bq = 3 ■ 10^^ T), and for other parameters as for 
the line 3 from Figs. 4 and 5, after solving Eq. (TT^ again, in Fig. 6 we present the wave 
frequency and the growth rate in terms of the temperature. Observe that, for example, 
at the normalized temperature T = 0.4 (the point C) the plasma-/? = O.Sme/mi, so due 
to the stronger magnetic field the mode is now in the electrostatic regime (compare with 
the point B from Fig. 5). At the same point we have the frequency and the growth rate 
u = 13.3 + iO.22 Hz. Hence, for such a stronger magnetic field the growth rate becomes 
about 50 times lower, as compared to the electrostatic case for the point A from Fig. 5. 

The points A and C in the previous examples may belong to the same magnetic loop. 
However, different loops may have different plasma-/? and these points may also represent 
such a situation. Therefore, the heating in different loops may be with or without a 
detectable variation of the magnetic topology. Observations of strong energy release events 
in the past (even in the range of flares) have shown that both scenar i os are indeed possible 



examp l es without magnetic va riations can be seen in 



mm; 



Pudovkin et al. 



JanssensI (1l972l ): 



Mayfleld fc Chapman 



(119981 ). The qualitative analysis described above will additionally 
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be supported below by some more numbers. 



According to iMcChesney et al.l ( 119871 . Il99ll ): ISanders et al.l (Il998l ) the stochastic 
heating by the drift wave is in action provided a strong enough wave potential amplitude, 
more precisely if 



mik^cj) 



2 ^2 _ _e0i 

't^iO 






> 1, 



and the maximum achieved ion velocity due to this heating is given by 



[eypie(f)i/{KTi) + 1.9]ni/ky 



(15) 



(16) 



The effective stochastic temperature is then T^ax = ''^i'^max/(3'^)- From Eq. ( 1T5|) it 
follows that the condition for the stochastic heating will sooner be satisfied in regions of a 
weaker magnet ic field, regardless of the start in g temperature . The physics of the heating 



is described in 



McChesney et al 



(119871. 



application to the solar corona in 



1991 



Ik 



Sanders et a 



Vranies fc PoedtsI (l2009a 



998l ) and its electrostatic 



d|), so the details of this 



will not be repeated here. We shall only stress that the heating is essentially due to the 
ion polarization drift. We shall apply these expressions using the ES and EM growth rates 
given above in order to quantitatively verify the differences due to eventual magnetic nature 
of the perturbations. 

For the parameters corresponding to the point A from Fig. 5, the expression ( 1T5|1 for 
a = 1 yields the required potential (pi = 61 V. The maximum achieved stochastic velocity 
from Eq. (ITB]) is 221 km/s, and the achieved stochastic temperature is T^ax = 1-97 ■ 10® K. 
Assuming some small accidental initial perturbations with the amplitude e(j)/{K,Ti) = 0.01, 
i.e., = 0.86 V we can calculate the time tg required to achieve the required value for 
the heating 0i = (f)exp{jtg). This yields tg = ln(0i/0)/7 = 0.4 s. Note also that here 
e0i/(/tTj) ~ 0.7. The total released energy density is Emax = ''^o^^j'^max/^ = 0.04 J/m^, 
and the energy release rate T^ax = Emax/tg = 0.1 J/(m^s). Hence, T^ax is about 1700 
times the required value for the coronal active regions [that amounts to ~ 6 ■ 10^^ J/(m'^s) 
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( iNarain fc UlmschneiderlllQQOl )] 



Taking as another example the point B in Fig. 5 (i.e., the same magnetic field 10^^ T 
but different density and starting temperature) yields Emax = 0.4 J/m.^, tg = 1.3 s, and 
^max = 0.3 J/(m^s), Tjnax = 1-97 ■ 10^ K] f3 = 2.5me/mj. So the present case is weakly 
electromagnetic and it is accompanied with the increase in the energy density and the 
energy release rate (because the density is higher), as compared to the electrostatic case 
from the point A, although it implies a longer growth time. The achieved stochastic 
temperature Tmax and velocity Vmax are the same as in the point A because a is kept fixed, 
and 01 = 61 V , = 0.34 V as in the previous case. 



On the other hand, taking as example the point C from Fig. 6, and the threshold (TTSjl . 
yields: = 0.34 V, 0i = 546 V, v^ax = 663 km/s, T^ax = 1-8 • 10^ K, E^^x = 3.67 J/m^, 
tg = 33 s, and T^ax = 0.11 J/(m^s). The stronger necessary potential here is due to the 
increased value of the magnetic field, see Eq. f fTSj) . Hence, the energy release rate Tmax is 
almost the same as for the point A, yet the characteristic time tg for the point A is more 
than 80 times shorter. In the same time, the maximum released energy density in the area 
with such a stronger magnetic field Bq is for about a factor 90 larger in comparison with the 
point A, with the achieved stochastic temperature that goes to 18 million K. The reason for 
the larger energy density is clearly the larger maximum stochastic velocity in the area where 
both the magnetic field and density are larger. The fact that plasma-/3 for this stronger 
magnetic field is only 0.3r?7,e/mj tells us that the increased stochastic energy density can be 
related to the EM effects and the coupling with the Alfven wave. However, a much more 
pronounced effect on the heating should be attributed to the increased magnetization of the 
plasma species. Thus, the areas with stronger background magnetic fields are subject to 
stronger stochastic heating. The magnetic field used here is in agreement with observations 
of active regions showing the field strength of a few times 0.01 T, that in fact may easily 
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go above 0.1 T f lLee et al.lll998l : ISolankill2003[ ). implying a possibly still stronger heating 



within the scenario presented above. 

Note also that the two potentials for the points A and C, 61 V and 546 V, respectively, 
are obtained assuming a = 1 in Eq. (!T5|) . These two potentials yield the electric field 
ky(j)i in the perpendicular direction 0.77 kV/m and 6.9 kV/m, respectively. Now, to have 
the threshold a = 1, the required potential 0i ~ \Bq, so a slight increase in these two 
parameters will yield even stronger electric fields. Taking as example Aj^ = 2 m, i?o = 4- 10^^ 
(instead of A^ = 0.5 m, Bq = ?> ■ 10^^ as in the point C) yields the perpendicular electric 
field at which the stochastic heating takes place Ey ^ 27 kV/ni. The three obtained 
values for the electric field, together with the corresponding magnetic field values, yield 
the E X S-drift (= E/ Bq) of the plasma as a whole in the perpendicular direction 77, 230, 
and 675 km/s, respectively. In view of the meter-sized perpendicular wavelengths these 
plasma fiows (drifts) could eventually be observed only by spectral analysis. Hence, we 
conclude that i) exceptionally strong perpendicular electric fields are expected during the 
proposed stochastic heating, and this particularly within stronger magnetic structures, and, 
ii) the perpendicular stochastic heating, as single particle interaction with the wave, is 
accompanied with collective plasma drifts. 



7. Summary and conclusions 

The results presented in this work could be summarized as follows. The kinetic theory 
of the drift wave shows that the mode is almost always unstable due to purely kinetic effects, 
and it couples naturally to the Alfven wave. The higher the plasma-/3 is, the better the 
coupling. The electrostatic drift wave in the solar corona is expected to be more unstable 
as compared to the regime in which the two modes are coupled. Essential for plasma 
heating is the electrostatic part of such an electromagnetic drift-Alfven wave. The heating 
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is stochastic by nature and, as shown in our previous works lVranjes &: Poedtd ( l2009al JblJdJdl). 
it possesses such properties that it is able to satisfy numerous heating requirements in the 
solar corona. From the analysis it also follows that the regions with stronger magnetic 
fields will be subject to much stronger heating. Note that such a relation between the 



temp erature and the magnetic field has been established long ago ( Ivan Speybroeck et al 



19701 ): in the present work we give a new, alternative explanation for this phenomenon. 



The electric field associated with the drift wave implies the possibility of the acceleration 
of plasma particles (primarily electrons) in the direction parallel to the magnetic field 
vector, and the development of drifts (in the perpendicular direction) of the plasma as a 
whole due to t he E x -B-drift tha t is the same for both electrons and protons. This issue 



is discussed in 



Vranies fc PoedtsI fl2009b 



jjd). The mean free path of the plasma species 
j is proportional to f^ • and therefore the parallel acceleration by the electric field is 
always more effective on particles that are already faster, i.e., those from the tail in the 
starting (Maxwellian) distribution, because those have more time/space to interact with 
the field. This will consequently result in a very different distribution function with a 
much longer high-velocity tail and resembling the ^-distribution observed in the outer solar 
atmosphere. The proposed electron acceleration within the present model appears as a 
natural development of the drift wave instability for which the source is clearly identified, 
thus removing the standard problem of various acceleration schemes that typically suffer 
from a common problem, the lack of a proper source. 

Some phenomena that follow fro m the presented s t ochasti c heating are not discussed 



here, but they are given in detail in 



Vranjes fc PoedtsI ( l2009bl ). These include the fact 



that the proposed model explains the better heating of heavier particles (i.e., heavier 
ions are better heated than lighter ones, while the ions in general are better heated than 
electrons). This follows after analyzing the mass dependence of the stochastic temperature 
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Tmax introduced earlier, with the help of Eq. f lT6|) . Also the better heating in the direction 
perpendicular to the magnetic field vect or, and the associ a ted tem perature anisotropy 
T± > Til, is self-evident and explained in lVranjes fc Poedtj (l2009bf ) as a consequence of the 



polarization drift that acts primarily in the perpendicular direction and becomes important 
at perpendicular wavelengths close to the ion gyro radius. 

The value of e(j)i/{K,Ti) in general determines the importance of nonlinearities. It 
turns out that in the examples discussed in the text this quantity is not small so that 
the presented scenario, which follows from the linear theory, may change considerably. In 
addition, the effects of nonlinearity in the drift-wave theory are determined also by making 
the ratio of the nonlinear term (i.e ., the convective derivative in the momentum equation). 



and the leading order linear term (JHasegawa fc Sato 



1989(@l ). The result can be written as 



{kyLn){kyPg)[e(f)i/{KT)] = kyLna. Because L„ :^ A^,, the proposed stochastic heating will 
be accompanied by various nonlinear phenomena. The most important nonlinear effects 
expected here include nonlinear 3-wave interaction that implies the well known double 
cascade (transfer of energy of a large amplitude drift-wave towards both longer and shorter 
wavelengths), and the anomalous transport caused by drift wave turbulence. These effects, 
however important, require numerical simulations and will be studied elsewhere. 
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Fig. 1. — The drift wave frequency for several different values of plasma fi in terms of the 
perpendicular wave-length. 



-25- 



0.6 r 



0.4 - 



3 



0.2 - 



0.0 



|3=0 

p =0.6 rrig/mj 

p =13 rrig/mj 

p =32 nrig/mj 




0.2 



0.4 0.6 0.8 

Xy [m] 



1.0 



Fig. 2. — The normalized growth rates corresponding to the frequencies from Fig. 1. 
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Fig. 3. — The normalized growth rate of the drift-wave obtained from Eq. f lT2|) in terms 
of the parallel wavelength and for three plasma number densities hq = 10^^ m~'^ (full line), 
no = 6 ■ 10^^ m~^ (dashed line), no = 10^^ m"'^ (dashed-dotted line). 
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Fig. 4. — The drift wave frequency obtained from Eq. flT2l) in terms of the plasma tempera- 
ture and for several values of the plasma number density (per cubic meter) . 
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Fig. 5. — The drift wave growth rate corresponding to the frequency from Fig. 4. 
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Fig. 6. — The drift wave frequency and the growth rate for i?o = 3 ■ 10 ^ T; other parameters 
are the same as for hne 3 from Figs. 4, 5. 



